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ABSTRACT 



The hypothesis that short GRBs arise from the coalescence of binary com- 
pact stars has recently gained support. With this comes the expectation that 
the afterglow should bear the characteristic signature of a tenuous intergalactic 
medium (IGM). However, fits to the observational data suggest that some de- 
tected afterglows arise in relatively dense gaseous environments rather than in 
the low density IGM. Here we show that considering the effect of red giant winds 
in the core of a star cluster may resolve this paradox if short GRB progenitors are 
contained in such an environment and close encounters rather than pure gravi- 
tational wave emission brings the compact objects together. Clear confirmation 
is provided here of the important notion that the morphology and visibility of 
short gamma-ray burst remnants are determined largely by the state of the gas 
in the cluster's core. 

Subject headings: stars: neutron — shock waves — globular clusters: general - 
gamma rays: bursts — hydrodynamics — stars: winds 



1. Introduction 



Although they were discovered roughly forty years ago (Klebesadel, Strong & Olson 



1973), Short 7- Ray Bursts (SGRBs) are still a mystery. All that can confidently be said is 
that they involve compact objects and highly relativistic dynamics. The most widely favored 
and conventional possibility is that they are produced by the coalescence of compact object 
binaries involving neutron stars and/or black holes, although other alternatives should also 
be considered (see e.g. Nakar 2007; Lee & Ramirez- Ruiz 2007, for reviews). Such systems, 



which are known to exist (Hulse & Taylor 1975 Burgay et al. 2003), are driven to coalesce 
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by energy and angular momentum losses to gravitational radiation. In this scenario, the 
compact binary would take hundreds of millions of years to spiral together, and could by 
then (especially if given a kick velocity on formation) have moved many kilo-parsecs from 



the site of its birth (Fryer et al. 1999 Bloom et al. 1999). However, fits to the observational 



data suggest that some detected afterglows (Berger 2007 Nakar 2007 Panaitescu 2006 



Nysewander et al.|2008[ ) arise in relatively dense gaseous environments rather than in the low 
density Intergalactic Medium (IGM). We must therefore remain aware of other possibilities. 
It may be wrong, for instance, to suppose that the binary is formed in isolation, since it 
could instead be the result of an encounter in a dense stellar system. 

Most stars in the Universe never interact strongly with others, at least during their 
adult life, after they have left the interstellar gas cloud in which they were born. However, 
there are various dense stellar systems such as star clusters and galactic centers, where stars 
are sufficiently close to their neighbors to make encounters significantly more likely. To get 
a sense of how crowded such a region is, consider the contents of the inner core of a Globular 
Cluster (GC). While in the Solar neighborhood the typical distance between individual stars 
is more than 1 pc, a sphere with a radius of 0.1 pc around the core of a GC contains over 
a hundred thousand, if not millions, of solar masses in stars following tight orbits about 
the nucleus. The stellar density is therefore more than a million times higher than that in 
the neighborhood of our Sun. In such environments, it is unavoidable that many single and 
binary stars will undergo close encounters and even physical collisions within their lifetimes. 
Indeed, it is thought that this type of interaction can lead to the formation of blue stragglers 



(Glebbeek, Pols & Hurley 2008), and possibly intermediate mass black holes (Portegies 



Zwart et al. 2004). It is in these environments, called dense stellar systems, that a compact 



object binary can be formed by two- and three-body encounters (Grindlay, Portegies-Zwart & 



McMillan||2006[ ). The subsequent merger, or possibly direct collisions of compact objects(|Lee 
& Ramirez-RHi|[2007l ), can produce a SGRB, and the resulting afterglow could then at least 
in part be due to the interaction of the relativistic ejecta with the stellar winds of the cluster 
members. Due to the large stellar density in the star cluster core, the interaction of the 
external shock can take place with a denser external medium than that of the IGM. Much 
of the effort herein will be dedicated to determining the state of the circumburst material in 
the cores of GCs, and describe how this external matter can affect the observable burst and 
afterglow characteristics. 
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2. SGRBs from Dense Stellar Systems 

First, we assume that different types of stellar objects are distributed homogeneously in 
a spherical core of radius r c = 0.1r C| _i pc, within which the number density, n c , and stellar 
velocity dispersion, a c , are constant. The gas density in the cores of GCs depends on the 
various types of stellar members that populate their interiors. In a typical cluster, winds 
from red giant stars tend to have the greatest effect on the core's gas density despite being 
significantly outnumbered by other stellar types. These winds are typically slow-moving and 
dense, with velocities on the order of v w = 10t> W) i km s _1 and mass loss rates between 1CT 7 
and 10 -6 M yr _1 . In a steady, spherically symmetric solution, the electron density is 

n w (r) « 3 x lOV^^Mw.-r^ 1 cm" 3 , (1) 

where /i e ~ 2 in a Helium gas, r = 10 16 r 16 cm and M w = 10 M W) _7 M yr _1 . The inner 
core of a GC contains iV* = 10 2 N Jft2 red giants, and so their mean separation is 

r± = 6.4 x 10 16 iV~ 2 1/3 r c _i cm. (2) 

We can thus estimate a minimum average density n± in the cores of GCs in the absence of 
gas retention as: 

n± ~ n w (r ± ) = SON^r^v'^M^rfi; 1 cm" 3 . (3) 

In this case, the SGRB would expand into a medium that is significantly denser than the 
IGM. 

For this discussion we will assume that the blast wave is adiabatic, i.e. its energy 
is constant with time, and effectively spherically distributed. This means that the energy 
{E = 10 50 i?5o erg) is the isotropic equivalent energy as, for example, derived from the 7-ray 
output. Deceleration due to the combined stellar winds starts in earnest when about half 
the initial energy has been transferred to the shocked matter, i.e. when it has swept up 
r _1 times its own rest mass. The typical mass where this happens is M& ec = E/(T 2 c 2 ) ~ 
5 x 10~ 8 E^iT^ 2 M Q . In the GC core, the mass within radius r is ^-p±r 3 , which gives the 
blast wave deceleration radius 

r d = 10^r^( T ^ Fi )- 1,3 c„, (4, 

A blast wave in a GC core thus decelerates at a much smaller radius than it would in the 
IGM. 



Under the assumption that energy conversion takes place primarily within the forward 
shock, the energy equation reads E hw (t) = ipT 2 M hw c 2 , where E bw is the blast wave's isotropic 
equivalent energy, ip is a constant of order unity ( Blandford fc McKee|1976 ), and M bw (r bw ) = 
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Ati J Q rbw n w {r)r 2 dr is the cumulative swept-up mass. The observed time interval during which 



most of the photons emitted (at a radius r bw ) are received is of the order of r bw / (AT c) (Sari, 



Piran & Narayan 1998). 



As usual, we assume that the dominant radiation process is synchrotron emission. If the 
energy in the magnetic field and electrons are taken to be a fraction e-Q and e e , respectively, of 
the thermal energy density, then approximating the electron energy distribution as a power 
law with an index p, the afterglow flux during the adiabatic expansion of the blast wave is 
given by 

Fy oc ^^/V^/X^M^e^M^ (5) 



for v sy < v < u c , and 



F v oc u'^-^E^M^-ht^e^ (6) 



for v c < v. The above relations are valid for varying energy or density, as can be seen by 
considering the implicit time dependence through E hw (t), n(t) and M hw (t). 

When the dominant variations are in the circumburst matter rather than energy, one 
may write t = (cj 1 AE)(M\ PN r\ )V] + J^ bw r 2 n(r)dr). Under this assumption, equations ([5| and 
(|gJ) reduce to F u oc M b L ~ p n (j3+1 ) //4 and F v oc M^~ p t~ l n( p ~ 2 ^ 4 , respectively. As the blast wave 
expands into the cluster winds, strong temporal variations compared to the canonical power- 
law decay can thus be produced as a result of changes on the properties of the star cluster, 
especially the mass-loss rates of the stars and their number density. The characteristic mass 
scale where this takes place is M bw > Md ec . This phase ends when so much mass shares the 
energy that the (3T < 1, setting a non-relativistic mass scale Mnr = E/c 2 ~ 5 x 10~ 5 E5qM & . 
Beyond this point, the event slowly changes into a classical Sedov- Taylor supernova remnant 
evolution. 

In the absence of characteristic scales in stellar ejecta and in the ambient medium, 
self-similar, spherically symmetric solutions exist, and they are widely used to interpret ob- 
servational data on afterglows. By contrast, the interaction of a GRB with a non- uniform 
medium is poorly understood. The presence of a density gradient will only affect the dy- 
namics of the GRB when the remnant size is comparable to, or exceeds, the scale length of 
the gradient. Before this time, the density can be treated as approximately uniform. Thus, 
it is only when > Tj_ that the details of the transition are important for the dynamics. 
This requires 

> 10 3 . (7) 



A close examination of equation ([7]) shows that if the winds of red giants are especially weak 
(i.e. M_7 < 1) or the burst's energy content is large (E 51 > 1), r± falls within the range of 
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relativistic expansion. Otherwise, this radius is sufficiently large that the interaction with 
the free wind of the nearest red giant is expected over the typical period of observation of 
afterglows. 

Depending on the wind properties of the stellar members as well as their number density, 
however, the density structure in this region could be quite complicated as the winds of 
stars interact at a typical distance ~ r±/2. Each free expanding wind encounters an inward 
facing shock, where the typical preshock densities are ~ An± ~ 320M w - 7 N^2 r c-i v w\ cm -3 . 
Kinetic energy is deposited in the shocked wind region in the form of heat, with temperature 
Tshock = (3m p /ie/16/c)(Av w ) 2 = 4 x 10 3 (At> w /10 km s -1 ) 2 K, where Ai> w is the speed of the 
material relative to the approaching shock. 

The shock interactions between nearby stars will be radiative if the cooling distance 
r coo i satisfies the condition k = r coo i/r± < 1. Using the cooling distances obtained by 



Hartigan, Raymond & Hartmann (1987) for shocks in the low-density regime derived from 
self-consistent pre- ionization, plane-parallel, steady models, we obtain k C 1. A cluster 
wind in this high cooling regime will have dense, cool structures between nearby stars and 
the SGRB blast wave could interact with them while still relativistic provided that < r±. 
Large-scale density inhomogeneities in the circumburst medium are, however, likely to result 
in distortions of the ejecta that might be more readily observed as the blast wave evolves 
into the non-relativistic phase. 



Cluster Winds Interaction Models for SGRBs 



3.1. Initial Model 



To build the initial model, iV* stars were randomly distributed in a three-dimensional 

1 /3 

volume such that the mean separation between stars is r_i_ ~ Av r c . The computational box 
is 8.0 x 10 17 cm in each dimension with the stars located within the central r r ~ 6.0 x 10 17 



cm. Based on the data from M15 (Dull et al. 1997 van den Bosch et al. 2006), the number 



density of red giants within r c is somewhere between 10 and 100, and for this particular 
simulation we set iV* = 10 2 , such that r± ~ 6 x 10 16 cm. Each star was given the same 
effective mass loss rate Mrg = 10~ 7 M yr -1 and stellar wind speed vrg — 10 km s" 1 . Based 
on Mrg and vrg, a spherical 1/r 2 wind profile was implemented for each cluster member 
(Figure [T]). The individual density profiles were then extrapolated so that they are in ram 
pressure equilibrium with the closest neighbor. 
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3.2. Dynamics of SGRB Remnants 

In the case depicted in Figure [T] the stellar winds of the individual members are dense 
enough to slow down the ejecta (for reasonable SGRB properties) to non-relativistic speeds 
before reaching r_|_, so that we expect the blast wave evolution as we see it to take place 
in the free-expansion phase. This encourages us to present a detailed account of resulting 
dynamics of the remnants within a dense stellar system following the onset of the non- 
relativistic phase. Common to all calculations is the initiation of the SGRB explosion as two 
identical blobs expanding in opposite directions into the circumburst medium. Calculations 
were done in three dimensions using the PPM adaptive mesh refinement code FLASH (ver 
2.5). The blobs and the circumburst medium are modeled as a cold ideal gas with 7 = 5/3. 
The numerical domain is a unprolonged cylinder in which the ejecta moves along the y-axis. 
In the inner region of each of the pancakes, the ejecta mass, M- s is distributed uniformly and 
for all runs we have used Tj = 2, Arj/rj ~ 0.6, 9j ~ 0.5, and rj ~ r d(rj). 

Without a detailed understanding of the exact shape and energy distribution of the 
ejecta, we have only an approximate description of how to construct the initial conditions. 



However, as clearly illustrated by Ayal & Piran (2001) and Ramirez- Ruiz & MacFadyen 



(2008), the late time evolution of the ejecta is rather insensitive to uncertainties in the 
initial conditions. We have considered various initial densities, angular widths, and shapes 
of the collimated ejecta and found that these are indeed unimportant in determining the 
late morphology of the remnant. This stems from the fact that at late times the mass of the 
remnant is dominated by the circumburst gas, which washes out any variations in the initial 
conditions of the ejecta. 

Detailed hydrodynamic simulations of the evolution of a GRB remnant in the GC core 
medium are presented in Figure [2j where the pressure contours of the expanding ejecta are 
plotted. As the gas collides with the external medium, a blast wave forms that propagates in 
the radial direction of motion and, over time, wraps around the various density discontinuities 
produced by the red giant stars. As it can be seen in panel (a) in Figure [2j the evolution 
of the blast wave suffers only minor distortions as the remnant sweeps a mass ~ E/c 2 . By 
the time the blast wave reaches the edge of the cluster's core, it becomes fairly distorted, 
though it still holds a jet-like appearance. Pronounced dimples are clearly visible in the 
lower portion of the blast wave showing where stars have been overtaken. Also shown in 
Figure [2] for comparison is the evolution of a spherical blast wave of comparable isotropic 
equivalent energy. 

The calculations above demonstrate how the dynamical evolution of SGRB remnants in 
dense stellar systems depends sensitively on the stellar wind properties of the cluster members 
as well as their density distribution. This confirms the notion that that the morphology 
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and visibility of short gamma-ray burst remnants are determined largely by the state of 
the gas in the cluster's core. The resulting evolution depends also fairly strongly on the 
properties of the GRB ejecta, especially its energy content, since it sets the non-relativistic 
mass scale for a particular stellar cluster. Figure [3] illustrates the dependence on energy 
of the remnant's morphology, where structures similar to those described in Figure [2] are 
clearly seen. A bow shock forms as each blob collides with the surrounding medium, which 
eventually wraps around the ejecta before the two expanding shells collide to form a single 
structure. However, in the low energy content case (panel a) the GRB remnant decelerates 
much more rapidly and although initially it may be highly nonspherical, the aspect ratio 
approaches unity as the two blobs expand and merge before reaching the edge of the GC 
core at about 2.3 x 10 8 seconds. Beyond this point, the evolution will follow the classical 
Sedov- Taylor supernova remnant solution. With higher energy content, on the other hand, 
deceleration to non-relativistic velocities occurs at much larger radii (panel c), and the blast 
wave is still moving through the core region before the two expanding shells collide to form 
a single structure. 



4. Discussion 

It is evident from the above discussion that the environment in the cores of GCs is a 
very rich one in terms of observable consequences. Even in the simplest case of red giant 
stars with identical wind properties, complex behavior with multiple possible transitions in 
the observable part of the GRB remnant's lifetime may be seen. Detailed, high resolution 
3D simulations of a spherical GRB exploding within the core of a dense stellar environment 
have been presented here, which show the resulting dynamics of the remnants following the 
onset of the non-relativistic phase. The resulting afterglow light curve will depend fairly 
strongly on the properties of the system, especially the mass-loss rates of the stars and their 
number density with various implications. On the one hand, it implies that one cannot be too 
specific about the times at which we expect to see transitions in the observed emission. More 
constructively, if and when we do see these transitions, they can be fairly constraining on the 
properties of the birth sites. If we continue to see the population of afterglows dominated by 
high inferred densities when compared to the IGM, this is support for blast waves in dense 
stellar systems, i.e. support for the origin of SGRBs from dynamically formed compact 
binaries or collisions. 

Stars in a dense stellar system interact with each other, both through their ionizing 
radiation and through mass, momentum and energy transfer in their winds. Mass loss 
leads to recycling of matter into the cluster wind gas, often with chemical enrichment. The 
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task of finding useful progenitor diagnostics is simplified if the pre-burst evolution leads 
to a significantly enhanced gas density in the immediate neighborhood of the burst. The 
detection of absorption signatures associated with the SGRB environment would provide 
important clues about the triggering mechanism and the progenitor. 
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Fig. 1. — Simulated red giant stars within the GC core. Each star is surrounded by a 1/r 2 
density profile, and is in ram pressure equilibrium with neighboring stars. The yellow surface 
is a contour of constant density (in M & c~ 3 s~ 3 units) around the star, and is not representative 
of the actual (unresolved) stellar radius. The size of the computational domain is (0.2 pc) 3 . 
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Fig. 2. — Morphology of the SGRB remnant 4.8 x 10 7 (a, top row) and 1.9 x 10 8 (b, bottom 
row) seconds after explosion. In the left panels, the pressure isosurface (in units of M c _1 s -3 ) 
is located at the outer boundary of the blast wave, and yellow/green/blue shading shows 
the density profile around the stars (in M Q c~ 3 s~ 3 units). In the early stages, the blast 
wave remains largely undisturbed. The size of the computational domain is (0.2 pc) 3 , and 
calculations were carried out in cartesian coordinates in three dimensions with six levels of 
refinement. The pressure slices in the right column show the region interior to the blast wave 
for each case. Inset panel: Behavior of a spherical explosion with a comparable isotropic 
equivalent energy = 5 x 10 48 erg. The initial high pressure region is in this case is a 
sphere in which the ejecta mass, and the thermal energy are initially uniformly distributed. 



Fig. 3. — The evolution of a SGRB remnant with varying energy at t — 2.3 x 10 8 , 10 8 and 
4.4 x 10 7 seconds (panels a,b and c, respectively). Shown are 3D pressure isosurfaces in 
normalized units (M c _1 s~ 3 ) for calculations with Eq = 5 x 10 47 , 5 x 10 48 and 5 x 10 49 erg 
(panels a,b and c, respectively). The blast wave has reached the edge of the core radius, at 
which point it has become extremely distorted. Pronounced dimples are clearly visible in 
the lower portion of the blast wave showing where stars have been overtaken. 



